Theoretical Reproduction of the Lineshapes of Nonlinear Conductance Observed in a 

Quantum Point Contact 
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The nonlinear conductance observed in a quantum point contact is theoretically reproduced for 
the entire range of applied bias. The single-impurity Anderson model with two reservoirs at different 
chemical potentials is studied for a sequential change of the gate voltage. The imbalance in the left 
and right Kondo coupling strength is introduced by the displacement of the Kondo impurity by the 
electric field produced in the constriction of quantum point contact. We reveal the origin of the side 
peaks and study the behavior of the height and width of the zero-bias anomaly. 



Understanding the nonlinear conductance in a system 
with nanocontact is indispensable in the study of na- 
noelectronics or nanomagnetism. The Kondo effect is 
often involved in studying its transport. A typical ex- 
ample of nanocontact is quantum point contact (QPC) 
or quantum wire [J-Q hi which an interesting sequential 
change of the dl/dV, where / and V denote the current 
and the bias, respectively, lineshape is observed with an 
increase in the gate voltage Vg- However, theoretical 
understanding for the observed nonlinear conductance is 
insufficient and unsatisfactory. The purpose of this study 
is to present the theoretical dl /dV lineshapes that show a 
similar sequential change of the dl/dV with Vg to the one 
reported in Ref. [3]. Our theoretical lineshape remark- 
ably fits the experimental data for the entire range of the 
applied bias used in the experiments. From this analysis, 
one may understand the physics of electron transport in 
a mesoscopic system in which a zero-bias anomaly (ZBA) 
and side peaks appear in the lineshape of the nonlinear 
conductance. 

In order to explain the observed nonlinear conductance 
in a QPC, the transport mechanism in a two-reservoir 
system with a mediating Kondo impurity must be inves- 
tigated when the system is under steady-state nonequi- 
librium (SSN). The issues of this study are two- fold: two- 
reservoir and SSN. The former duplicates all basis vec- 
tors describing the back and forth movements of the elec- 
tron. For this reason, additional resonant tunneling lev- 
els, which are the origin of the side peaks, occur. The 
latter shifts the Fermi level of the reservoir and allows 
only unidirectional resonant tunneling if the system is in 
the ground state. This unidirectional resonant tunneling 
makes the state of SSN different from that of equilibrium. 

We obtain the formula for dl/dV by differentiating 
the current formula provided in Refs. [f| and @ with re- 
spect to the bias. We fix the chemical potential of the 
Kondo impurity at the Fermi level of one of the metal- 
lic reservoirs. Then, the differential conductance at zero 



temperature is given by 
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where e, a, and fi denote electronic charge, spin, and 
the Planck constant divided by 27r, respectively. Here, 
Pi a (eV) is the local density of states of the Kondo im- 
purity at bias V. The factor T tot is given by the com- 
bination of featureless coupling functions, i.e., T tot = 
r L r R /(T L + T R ) 0, where the superscripts L and R 
stand for left and right reservoirs of a QPC, respectively. 
We omit the term containing dpi a (u>)/dV because its 
contribution to dl/dV is negligible. We show its neg- 
ligence in a different study. 

The spectral function of an up-spin electron, for ex- 
ample, at the Kondo impurity is given by the real part 
of the inverse of the matrix M r ||, i.e., Ptf(u) = 
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note the energy level of the Kondo impurity, Coulomb 
interaction, and the average number of down-spin elec- 
trons occupying the level respectively. All the matrix 
elements, except the eight ^/-elements, have additional 
self-energy terms, /3 m „[zS^(w) +iE R (u)] = 2(3 mn A, for a 
fiat wide band, where A = (T L + T R )/A and T R cx \V£ | 2 , 
where V R denotes the hybridization strength between the 
Kondo impurity and the right reservoir [7|. The coeffi- 
cient p m n is given in [9]. We use A to indicate the unit 
of energy. 

The matrix elements represented by 7 are given by Q 
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where the operator e^) denotes the fermion annihi- 
lation at the state-fc(i) in the lead (impurity), j77 = 

V^cl^Cii), and 5 means the deviation from the average. 
The element Uj T , on the other hand, is given by 



TT L - 
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One may roughly understand the roles of matrix ele- 
ments as follows: Jll(rr) represents the degree of Kondo 
coupling; "fLR and 7j represent connecting mechanism; 
The real part of [/-elements determines the position of 
Coulomb peak, whereas the imaginary part gives rise to 
the asymmetry in the dl/dV lineshape. A symmetric 
lineshape is obtained when (rii a ) = 1/2. 

The spectral function Pif(u>) given by the reduced ma- 
trix of Eq. ([2} has five peaks in general. The posi- 
tions of five peaks can be obtained from the zeros of 
the determinant of M r in the atomic limit: a Kondo 
peak at uj 1 = 0, two Coulomb peaks at u/ « ±C//2, 
and two additional resonant tunneling level peaks at 
u> = ±[( 7 £ L + jl R )/2 + ( 7LR 7J ) 2 + 0(U-i)}^ for 
a large U. One can also find the spectral weight of each 
peak. It is interesting to note that the additional tun- 
neling level can have a meaningful amount of spectral 
weight only when a bias is applied (7,/ 7^ 0). Therefore, 
we call the peak at the additional tunneling level the 
steady-state resonant tunneling level (SSRTL) peak. 

We have discussed the effect of bias, which allows only 
unidirectional resonant tunneling at the ground state, 
above. In Fig. 1, we present the graphical illustrations of 
jlr and 7j on the basis of the corresponding operator 
expressions, which are third order processes of hybridiza- 
tion. Both exchange (black or red @) and singlet hopping 
(blue or red @) processes in equilibrium are depicted in 
Fig. 1. Applied bias makes the second part of Fig. 1 neg- 
ligible when we calculate the expectation value. Thus, 
the condition Jlr ^ 7,/ is obtained. Note that only sin- 
glet hopping process can establish current. Therefore, 
the unidirectional resonant tunneling (without black @ 
on the left part of Fig. 1) dominates the Kondo process 
under bias. In this study, we do not calculate the values 
of 7. We will use them as free parameters within the 
condition 7^ > 7,/. 

In order to perform the calculation for the spectral 
function, we must know what happens in the constric- 
tion of a QPC when Vq changes. Calculations by spin 
density functional theory 13, 11 1 show that a magnetic 



FIG. 1: Motions of electrons indicated in ^lr (+ sign) and 
7.7 (— sign). Green sold dot indicates a Kondo impurity. The 
motion denoted by @ in red performs exchange process along 
with @ in black or singlet hopping with @ in blue. 



impurity is formed spontaneously at the center or a dis- 
placed position in the QPC constriction at a low Vq- 
When a bias is applied, an electric field will be created in 
the constriction and the magnetic impurity will be dis- 
placed to the direction opposite to the electric field, as 
shown in Fig. 2 (a). As a result, one has jll > Irr- 
As Vg increases, the magnetic impurity will move back 
to the center of the constriction because the number of 
electrons in the lead increases, as shown in Fig. 2 (b) 
{ill ^ Irr)- When the magnetic impurity reaches near 
the center of the constriction, we assume that jll = Irr 
is established, as shown in Fig. 2 (c). Upon reaching the 
center, the magnetic impurity stays there even at a high 
Vg- However, both jll and 7^ increase continuously 
as Vq increases further. 

According to the scenario discussed above, we deter- 
mine the specific manner in which the matrix elements 
(Ujf 1 and 7) and T tot change with Vg- One can see 

that the real part Re[J7j± fl ] is composed of the corre- 
lation in numerator and the fluctuation in denomina- 
tor, i.e., Re[J7j± fl ] = C^'±/F^'^ L . Since fluctuation is 
proportional to the number of particles N c , where N c 
is the number of conduction electrons in a QPC, and 

N c cx Vg [3, we obtain F^'± cx Vq 2 . We propose a 

L,R 
U± 



specific expression that shows the dependence of Fj 
on Vg- Fy±* = P\/Vg — q- The constants p and q are so 
determined that Fy±* = 1 for Vg = 1 and F^' ± = 6 for 
Vq = 17, as shown in the first two columns of Table I. 
The values of Vg indicate the steps in which Vg is varied. 

In the region where ^ll > Irr. (Figs. 2 (a) and 2 (b)), 
one must consider L — R unbalance in the strength of 



L R 

Kondo coupling. Only C n ' + that includes j 



■+L,R 
Hi 



would 

reflect the L — R unbalance. We consider linear changes 
in Vq for C v ' + up to Vq = 9. For Vq > 9, however, ^ll 
and "fRR are balanced and = Cff + = 1. We choose 
Vg — 9 as a border separating the L — R balanced and 
unbalanced region. We also consider linear relationships 
for Jll(rr) and T tot . We set ^lr. = 7./ for the unbal- 
anced region, i.e., Vg < 9. However, in the balanced 
region, Vq > 9, we consider nonvanishing exchange pro- 
cess in the second part of Fig. 1 and allow 7^ to increase 
further. 

We first fix the values of matrix elements of M r for 
Vg = 1 as shown in the lowermost row of Table I that 
afford the lineshape similar to the one obtained experi- 



3 



FIG. 2: Location of a magnetic impurity (black dot) under 
bias. The strength of Kondo coupling is represented by the 
line thickness, (a) Low Vg'- Jll > "/rr- Green dots denote 
electrons in the leads, (b) Vg before balance: ^ll ^ Jrr- (c) 
Vg after balance: 7ll = jrr- 



Table I: Change in parameters with Vg (7.7 = 0.5) 
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FIG. 3: (a) Sequential change in the dl /dV lineshape when 
Vg changes from 1 to 17. The curve in red (Vg = 9) is 
the border that corresponds to the 0.7 conductance anomaly, 
(b) Comparison of theoretical data (red line) with the ex- 
perimental data (crosses) of the closest-spaced side peaks in 
Fig. 1 (b) of Ref. [4]. Weak asymmetry is given by setting 
lm[Uj + ] = 0.05. We set /3n = 0.252, f3 12 = #21 = 0.246, 
P22 = 0.256, and other /3 m „ = 0.25. 



mentally for the lowest Vg- The sequential change of the 
dl/dV lineshape with Vg is obtained by using the val- 
ues in Table I that are determined by the Vc-dependence 
discussed above. We list the values only for odd numbers 
of Vg in Table I. Figure 3 (a) shows the characteristics 
of the sequential change of the nonlinear conductance in 
a QPC, such as the change of ZBA and the behavior of 
the side peak with an increase in Vg- One can observe 
two distinct groups of lineshapes, denoted by different 
colors. As Vg increases, the SSRTL peaks first move in- 
wards, and then move outwards. The behavior of SSRTL 
changes at Vg = 9, at which point the two SSRTL peaks 
are closest in bias. The values listed in Table I depend 
on top-gate voltage and sample characteristics. Different 
sample or top-gate voltage affords different table and the 
corresponding set of dl/dV lineshapes. 

A remarkable comparison with experiment is shown in 
Fig. 3 (b). We choose the experimental dl/dV curve that 
is the one with the closest-spaced side peaks of Fig 1 (b) 
of Ref. [4] for comparison. Our theoretical lineshape is 
obtained by setting 7ll = Irr = 0.5, 77 = j LR = 0.65, 
Re[Ujf] = 1.17, Im[U$ + ] = 0.05, and T tot = 0.62. Since 
the experimental curve is slightly asymmetric and the 
peak is a little bit offset from V = 0, we introduce a small 
nonvanishing Im[E/j + ]. We discuss about asymmetry of 
the lineshape more below. Perfect fitting is obtained, as 



FIG. 4: Characteristics of ZBA. (a) ZBA height (A/i Z ba) vs. 
G max . (b) Full width at half maximum (FWHM) of ZBA vs. 
Gmax- The dots are obtained from Fig. 3 (a) and the crosses 
are the data of Ref. [4] for Ko P = +4 V. The units of A/izba 
and G are 2e 2 /h. Inset: Conductance vs. gate voltage. 



shown in Fig. 3 (b). 

The ZBA peaks shown in Fig. 3 (a) provide us with 
another interesting comparison with the experiment. In 
Fig. 4, we present the height and width of ZBA vs. 
G max , the maximum value of the corresponding ZBA, 
along with the experimental data reported in Ref. [3| 
for Vt p = +4 V. Our theoretical data given in Fig. 3 
(a) reproduce the ZBA characteristics quite well. One 
may understand the peculiar behaviors shown in Fig. 4 
in terms of the competing relationship between the ZBA 
and the SSRTL peaks as Vg changes. In the region 
Vg < 9, increasing Vg has two effects, i.e., recovering 
the L — R symmetry of Kondo coupling strength and 
enhancing the Kondo peak by suppressing the fluctua- 
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tions ((Sj^ L ' R ) 2 ) 1 ^ 2 - These two effects are conflicting in 
determining the ZBA peak because the former causes in- 
ward movement of the SSRTL peaks with taking spectral 
weight from the ZBA. These conflicting roles of Vq lead 
to a local maximum of AhzBA in the region Vq < 9 in 
Fig. 4 (a). However, in the regime Vg > 9 in which the 
L — R symmetry is retained, only one role of Vg , enhanc- 
ing the Kondo peak, remains. Interestingly, however, the 
unitary limit of conductance prevents the height of ZBA 
from increasing further. As a result, another local maxi- 
mum in Fig. 4 (a) appears in the region Vg > 9. Abrupt 
increase in the width of ZBA from Vg — 9, as shown in 
Fig. 4 (b), is also resulted from this unitary limit. We also 
plot the maxima of Fig. 3 (a) vs. Vg in the inset Fig. 4 
(b). A slight curve bending is shown near 0.7(2e 2 //i). 

The dl / dV lineshape becomes asymmetric if the imag- 
inary parts of do not vanish. Since lm\U ^ ] oc 
(1/2 — {riiu)) as shown above, the asymmetry in the line- 
shape indicates the existence of spin polarization [13j in 
the magnetic impurity, i.e., (n^) ^ (nil). Therefore, the 
Kondo effect and spin polarization may not be incompat- 
ible. In order to identify the origin of spin polarization, 
a deeper analysis for the behavior of magnetic impurities 
created in the constriction of a QPC. Recently, Song and 
Ann [14| proposed a possibility of ferromagnetic coupling 
between the magnetic impurities. An extended Ander- 
son model considering the ferromagnetic coupling may 
be suitable for explaining both the Kondo and spin po- 
larization simultaneously. 

In conclusion, we obtained the dl/dV lineshapes ob- 
served in a QPC in terms of the Anderson model with 
two reservoirs. We identified the side peaks of the dl/dV 
as the SSRTL peaks. We showed that the sequential 
change of the dl/dV lineshape with Vg is separated into 
two parts (Figs. 3 (a)) and the border corresponds to 
the point of 0.7 anomaly (Inset of Fig. 4 (b)) at which 
the curve characteristics of ZBA height and width are 
apparently different (Figs. 4 (a), (b)). We also identi- 
fied that the strengths of the L — R Kondo coupling are 
just balanced at the point of 0.7 anomaly (Table I and 
Fig. 1 (c)). The theoretical skill used in this study can be 
extended to other mesoscopic systems with a mediating 
Kondo impurity between two reservoirs under bias. 
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